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Abstract 


ISEE-1 low energy electron measurements have revealed the frequent 
presence of field-aligned fluxes of few hundred eV electrons in the 
geomagnetic tail lobes. In the northern tail lobe these electrons are most 
prominent when the interplanetary magnetic field is directed away from the 
sun. This characteristic helps identify the electrons as polar rain electrons 
that are usually identified by their uniform precipitation pattern over the 
earth's polar caps. Although tail lobe electrons have been seen previously at 
high altitudes, the ISEE data provide the first demonstration of their field 
aligned nature and this characteristic places a stringent constraint on their 
origin. By mapping the tail lobe velocity distribution function into the 
solar wind, we have confirmed previous suggestions that the polar rain is 
indeed of solar wind origin and is due to the access of electrons to the 
magnetotail lobe. More specifically, however, it is demonstrated that the 
more energetic component of the polar rain is composed of electrons from the 
solar wind "strahl" - a field-aligned component of the solar wind which is 
difficult to measure but which is thought to be caused by the collisionless 
transit of hundred eV electrons from the inner solar corona to 1 AU. The 
preferential entry of the strahl into that polar cap most directly connected 
magnetically to the sun explains a known north-south asymmetry in polar rain 
Intensity. Furthermore, it is suggested that cases of higher energy polar 
rain that have been thought to require a magnetospheric field aligned 
acceleration mechanism are more likely due to an enhanced intensity of the 
solar wind strahl caused by easier access of solar corona electrons to 1 AU 
and the polar caps. This idea is supported by the observation that previously 
reported cases of intense polar rain are invariably associated with tenuous or 
high speed solar wind that is expected to produce enhanced strahls. It is 
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shown how polar rain measurements may provide a good measure of (1) the 
temperature and the state of statistical equilibrium in the solar corona, and 
(2) the degree to' which interplanetary field lines connect with those of the 
earth. 

1. Introduction 

Although magnetospheric field lines poleward of the auroral zone have 
long been known for their dearth of plasma, sensitive measurements in recent 
years have revealed the presence of low energy electrons on these field lines. 
These few hundred eV electrons frequently precipitate rather uniformly over 
the entire polar cap and have been termed the "polar rain" by Winningham and 
Heikkila [197^3. They exist on field lines that lead to the lobes of the 
magnetic tail and indeed measurements in the northern lobe [Yeager and Frank, 

i r 

1976] have revealed the presence of these electrons at =3-7R e . 

An interesting characteristic of these high latitude electrons is that 
their access to the magnetosphere is regulated by the polarity of the 
interplanetary magnetic field (IMF) sector structure. In the northern 
hemisphere of the magnetosphere, 300-500 eV electrons were found to be much 
more intense when the interplanetary field was pointing away from the sun 
(away sectors) [Yeager and Frank, 1976]. At low altitudes over the polar caps 
a north-south asymmetry was often detected with the north to south polar cap 
intensity ratios tending to be greater or less than unity depending on whether 
the earth was immersed in an away or toward sector [Fennell et al., 1975; Meng 
and Kroehl, 1977; Mizera and Fennell, 19783. It has also been noted that 
spectra of these electrons are similar to the cusp [Winningham and Heikkila, 


1974] and the magnetosheath [Yeager and Frank, 19763. Fennell et al, [1975] 
presented one example of rather good agreement among nearly simultaneous 
measurements in the tail lobe, over the polar cap, and in interplanetary 
space. 

The above facts led to the suggestion that low energy electrons have 
unimpeded access to the polar cap [Winningham and Heikkila, 1974; Fennell et 
al., 1975] in the same manner as energetic solar flare electrons [e.g, 
Paulikas, 1 974 ] . This concept is illustrated in Figure 1 which shows field 
lines from one polar cap connecting to interplanetary field lines leading 
directly back to the sun while field lines from the other hemisphere connect 
to field lines continuing far out into the solar system. Which polar cap 
connects most directly to the sun is controlled by the interplanetary sector 
structure. A directed flux of electrons moving along field lines and coming 
from the sun (open arrows) would preferentially be directed into one polar cap 
or the other in order to produce the north-south asymmetry. The 
interplanetary heat flux was suggested as the source of this asymmetry 
[Fennell et al., 1975], but detailed comparisons of interplanetary 
distributions functions and magnetospheric measurements have not been carried 
out. 

Other measurements have also shown that the intensities over the polar 
cap at any given time often exhibit large scale nonuniformities. Meng et al. 
1977 have presented evidence for a dawn-dusk intensity gradient that depends 
on the IMF. When the interplanetary is positive (negative) the northern 
hemisphere fluxes tend to be higher near dawn (dusk). The situation reverses 
in the southern hemisphere. These variations are similar to IMF dependence of 
the polar cap electric field [Heppner, 1972; Burke et al., 1979] with regions 
of strong E corresponding to regions of high precipitating flux. 
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Corresponding dawn-dusk variations of electron intensity were not seen in 
the tail [Yeager and Frank, 1976] but this may be due to the inability of a 
more distant 3paoecraft to sample different spatial regions in times that are 
short compared to the variability of polar rain. Very recently Cussenhoven et 
al. [1984] have discovered a day to night decrease (increase) in integral flux 
(average energy) in low altitude measurements of the polar rain. 

In their original study of ISIS traversals of the polar cap, Winningham 
and Heikkila [1974] defined two additional classes of somewhat more intense 
and energetic polar precipitation. "Polar showers" were narrow struotures 
embedded within broader polar rain regions which undoubtedly correspond to 
polar cap auroral arcs [see also Hardy, 1984], These structures will not be 
of interest here. "Polar squalls" were broader regions of enhanced keV 
electron precipitation associated with magnetically disturbed intervals. Low 
energy photoelectrons emitted from the polar cap were often seen returning to 
low altitudes within these squall regions whereas within polar rain regions 
these particles were usually moving outward and presumably escaping along open 
field lines [Winningham and Heikkila, 1974? Winningham and Gurgiolo, 1982], 
These returning photoelectrons led to the suggestion that parallel electrio 
fields over the polar cap were involved both in accelerating the more 
energetic squall electrons and causing photoelectrons to return to low 
altitudes. Heppner et al. [1981] have presented further evidence for such 
parallel electric fields. 

Foster and Burrows [1976] also studied intervals of uniform keV-type 
polar precipitation which were similar to polar squalls except that each of 
their two events aovered practically the whole polar cap and each occurred 
during the quiet interval following a 3torm. Meng and Kroehl [1977] presented 
two similar events, Foster and Burrows [1976; 1977] also suggested that these 


more energetic electrons were accelerated by field aligned electric fields. 
This proposal was based largely on the facts that (1) the north-south 
asymmetry eliminated the possibility that the electrons were bouncing on 
closed field lines, and (2) simultaneous interplanetary menaurements failed to 
reveal oomparable fluxes that might have been entering the magnetosphere, 

Meng and Kroehl had no interplanetary measurements and they favored a solar 
origin for their intense precipitation. Comparable measurements of intense 
fluxes in the tail lobes have not been reported. Currently, a field aligned 
acceleration mechanism seems to be favored a3 the explanation for these more 
energetic polar rain-like events seen over the polar caps [Foster and Burrows, 
1977 ; Mizera and Fennell, 1978 ], The recent measurements (Gussenhoven et al., 
1984 ] even suggest a role for the magnetosphere in modulating normal polar 
rain. 

In the following sections we describe ISEE-1 measurements of anisotropic, 
field-aligned, tail lobe electrons whose occurrence characteristics identify 
them as the tail lobe component of the electron distribution function that 
precipitates as polar rain. We then examine the question of solar wind origin 
of the polar rain U3ing contemporary knowledge of anisotropic solar wind 
distribution functions that was not available at the time of earlier 
magnetospheric measurements. By mapping these distribution functions along 
field lines between the solar wind and the tail lobe, we quantitatively 
demonstrate that the polar rain is indeed of solar wind origin with the 
north-south asymmetry due to the previously unknown anisotropic component of 
the solar wind - tile strahl. We review information on the solar origin of the 
strahl and reinvestigate the question of whether intense polar rain events 
could be of solar origin rather than due to magnetospheric acceleration. We 
conclude that the limited information on these intense events is consistent 
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with this solar wind hypothesis although the possibility of further 
magnetospheric acceleration remains a possibility. The solar origin of the 
polar rain has various implications for both solar and magnetospheric physics 
whioh are discussed. 


2. Instrument and Data Reduction 

The ISEE 1 electron spectrometer was designed primarily to provide high 
time resolution, three dimensional measurements in the relatively dense plasma 
at the magnetopause and beyond [Ogilvie et al., 1978], To achieve this goal 
six individual sensors with 8.5° x 11° viewing angles were mounted to form two 
mirror imaged triaxial spectrometers, each with three mutually orthogonal 
fields of view. The six 3en3ors are simultaneously sampled at logarithmically 
spaced energies extending over a range which is set to either 11 eV to 2.06 
keV (Mode 2) or 109 eV to 7.28 keV (Mode 3) when the spacecraft is in the 
magnetotail. Each energy sweep takes 0,5 s or 1/6 of a 3 second spacecraft 
spin. As the spacecraft spins the individual sensors, which make angles of 
±16,2, ±31.7i and ±53.5 degrees with the spacecraft equatorial plane, sweep 
out 6 spatial bands on the unit sphere while they are each undergoing 6 energy 
sweeps, Thu3 each energy is sampled 6 times by each of 6 detectors giving 36 
widely separated points on the unit sphere, half of which are located 
diametrically opposite to another simultaneously sampled point. In the high 
(low) bit rate mode every third (sixth) 3 3 spin is sampled giving complete 
data 3ets every 9 or 18s. These measurements yield 36 point determinations 
of the pitch angle variation on every spin when combined with the 1 minute § 
data from the experiment of C. T. Russell obtained via the ISEE data pool 
tapes. 
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In the course of calculating moments of the electron distribution 
function, f, the valuo of f i3 determined for each of the 36 measurement 
directions at each energy. To obtain a rough measure of the electron 
isotropy, at each energy we calculate an average f A for all measurements with 
pitch angles 90° ±30° and an average f„ for all the remaining f values having 
pitch angles within 60° of the parallel or anti-parallel field direction. 

Plots of f„ and f A versus time have proven useful in identifying magnetotail 
anisotropies. 

3. Polar Rain in the Tail Lobe 

Initial scans of ISEE 1 electron data in the magnetotail revealed the 
presence of occasional intervals of persistent anisotropic field aligned 
fluxes in the tail lobes with energies typically ranging up to a few hundred 
eV. An example of suah anisotropic fluxes is shown in Figure 2 where the 
measured log of the distribution function f (irregular solid lines) is plotted 
versus pitoh angle. Each tic mark on the ordinate represents one decade in f 
and the adjacent energies have been offset 1.25 decades to prevent 
overcrowding. Fluxes are low, as is characteristic of the tail lobes, so 
instrumental counts have been averaged over 30 consecutive spins (five 
minutes). Only background counts are detected for pitch angles near 90° for 
energies above 180 eV, but the counting rates rise to several tens of oounts 
per spin above background for the lowest and highest pitch angles. One count 
levels are indicated by dashed lines. Note that the logarithimic scale of f 
tends to deemphasize this more dramatic increase in the counting rate near the 
field direction. The fluxes are approximately symmetric parallel and 
anti-parallel to the field direction and they are only detected within about 
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30° of the field direction. Another important characteristic illustrated by 
Figure 2 is the energy dependence of the anisotropy; the lower energy 
electrons are invariably more isotropic than those at higher energies. 

To investigate the persistence of thc<se anisotropic electron fluxes we 
utilize plots of the average f„ and f A as defined above. Since f„ i 3 an 
average of all measurements within 60° of the field direction, it ha3 a 
magnitude that is substantially diminished relative to the actual f value near 
0° or 180° pitch angle. In spite of this fact, these average f„ and f A values 
serve as useful indicators of the anisotropy. An example of suah data 
oovering six hours on April 18, 1978 is shown in Figure 3. The lines 
representing f„ and dots representing f A are plotted on a logarithimic scale; 
every third point is omitted to prevent overcrowding and successive energies 
are offset by one deaade. The example of Figure 2 comes from this interval 
and the separation of the traces in Figure 3 can be seen at the energies where 
the anisotropy is apparent in Figure 2. The fluxes are relatively steady on 
thi3 scale of hours until the plasma sheet (more intense energetic electrons) 
i3 approached shortly after 0700. 

The observed characteristics of these lobe electrons suggest that they 
may be polar rain electrons. The fact that anisotropic lobe electrons are 
seen in the northern tail lobe on some days but not on other days is 
consistent with the preferred-polar cap entry process. Their steady nature 
when present is similar to low altitude observations and their energy range is 
similar to polar rain. Although previous measurements on tail lobe field 
lines did not reveal any pitch angle anisotropy [Yeager and Frank, 1976], 
there appear to be good reasons for this. These workers were using a wider 
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angle deteotor -hat summed counts over 55° of azimuthal spin. They wore also 
measuring primarily at 3-7 R Q where the anisotropy would be less pronounood 
than in the deeper tail (see discussion below). At these altitudes the field 
direction was also probably not often directed into their viewing angle. 

To strengthen the evidence that anisotropic lobe electrons are indeed 
polar rain electrons, a statistical survey was undertaker* of tail lobe hourly 
averages. Various quantities were averaged including f„, f A and the ratio of 
their hourly averages. Using the electron measurements along with the data 
pool magnetic field and occasionally the ion spectrograms of L, A. Frank, 399 
hours were identified when ISEE 1 was in the tail lobe (excluding magnetopause 
boundary layer and planma mantle) during the interval February 10 - June 7* 
1978. During these hours the lobe electrons were Investigated a3 a function 
of interplanetary sector structure or interplanetary B y component. Of these 
399 lobe hours, only 208 (52%) corresponded to hours when IMF measurements 
were available and only 43 of these corresponded to positive when lobe 
electrons are expected to be present. If, however, the sector structure as 
determined from ground magnetograms was used [Lincoln, 1983], 351 (88%) of the 
hours could be associated with toward or away sectors and 89 of these hours 
corresponded to away sectors. Although the direct measurement is obviously 
more accurate when available, a comparison of the two methods confirmed the 
relatively high accuracy of the ground determinations and we use the latter 
method because of its greater completeness. 

In Figure 4 we show a histogram of the parallel to perpendicular f ratios 
for the two sector polarities. It is important to realize that not only is f„ 
a lower bound on the aotual maximum f at 0° pitch angle due to averaging over 
a range of pitch angles, but is an upper bound set by the instrumental 
threshold which is undoubtedly much higher than the actual f at 90° pitch 
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angle. Both effects imply that the computed ratios are lower bounds. In 
toward sectors when low energy electrons are expeoted to be absent, the 
averages of f A and f„ are both near their instrumental background level and 
their ratio is greater than 1. 1 only 131 of the time. During away sectors 
where polar rain electrons are expected to be present, the ratio is greater 
than 1, 1 9036 of the time. The absolute average intensity for f„ was 
1, 65x1 0 -29 s^ cm" 4 * 6 for away sectors at 180 eV, a value 2.4 times higher than 
for toward sectors. This factor undoubtedly would have been higher except for 
thu fact that low values in toward sectors were limited by the instrumental 
threshold. 

Figure 4 is clearly in agreement with the previously reported polar rain 
dependence on interplanetary sector structure and we ocnclude that the 
anisotropic lobe electrons are indeed the magnetotail component of the polar 
rain. We prooeed to investigate the relationship of these electrons to the 
solar wind. 

4, Interplanetary Origin of the Polar Rain 

Although access of solar wind electrons ha3 been suggested as the source 

for the polar rain, little consideration has been given to whether asymmetries 
in the observed solar wind distribution functions are quantitatively 
compatible with the north-south polar cap asymmetry. Now, the additional 
observation that the magnetotail source of polar rain electrons is highly 
anisotropic puts a new, stringent constraint on the solar wind origin 
hypothesis. Rather than simply requiring that the average solar wind flux be 
consistent with magnetosphere observations [Fennell et al., [19753, or more 
specifically that the field aligned solar wind electrons are consistent, we 
mu3t now ask whether complete solar wind distribution functions, when mapped 
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along field lines into the higher fields of the magnetotail, are consistent 
with the observed anisotropy. A? will be illustrated below, such 
compatibility i3 far from evident since theory indicates that the particle 
distributions tend to beoome more isotropic when mapped to a stronger field 
region. 

4,1 Happing Interplanetary f 

The solar wind distribution function f can be represented as the sum of 
three maxwellian or bi-maxwellian components. 

f =f 0 *f h .f a to 

The "core", f c , and "halo", f^ , are well known components of the solar wind 
[Montgomery et al.,1968]. The core, i3 a cool component with temperature 
=1x10 5 °K which accounts for about 94* of the solar wind density [Feldman et 
al., 1975]. This cool component i3 important only at very low energies (<!40 
eV) that ars either below our measurement range or in channels contaminated by 
the 3pac<i,«'.raft plasma sheath; we ignore this component in the present 
discussion. The "halo" component, f^ , can be approximated a3 a maxwellian 
with higher temperature and lower density, 

i’h = V <^V) exp(-v 2 -v a 2 )/ u 2 (2) 

where u 2 = 2kT h /m e and T h = 8x1 0 5 . 

The "strahl", f , is a field aligned component which was not discovered 
until measurements from relatively narrow angle detectors on the Helios 
spacecraft were surveyed [Rosenbauer et al., 1976; 1977]. Earlier surveys had 
used detectors with apertures that were large compared with these field 
aligned electrons and hence the Intense directed fluxes were indistinguishable 
from weaker broadly distributed fluxes. The strahl has subsequently been 
confirmed by Feldman et al. [1978] and Ogilvie and Scudder [1981]. It is 
thought to be due to electrons of the solar corona which reach the vicinity of 


1 AU wit.i row collisions and minimal scattering [R03enbauor et al., 1976; 

1977; Olbert, 1981] (see section 6). These electrons are very field aligned 
because they tend to conset ve their first adiabatic invariant (u s v.'VB) as 

i j' * 

they move away from the sun in the decreasing solar (interplanetary) field. 

The lower energy strahl electrons are scattered more than those at higher 
energies due to the strongly energy dependent coulomb collision frequency, and 
therefore lower energy strahl electrons are observed to be more isotropic than 
those at higher energies [Rosenbauer et al., 1976; Olbert, 1981} 1983], The 
halo distribution may in fact be simply the lower energy corona electrons that 
have been scattered enough to make them fairly isotropic [Rosenbauer et al,, 
1976; Scudder and Olbert, 1979a, b], A more detailed discussion of the escape 
of strahl electrons from the sun will be given in section 6, For the present 
we note that these electrons are important in carrying the interplanetary heat 
flux and they have a preferred direction of motion away from the sun and along 
the local magnetic field. This motion gives them the potential for 
preferentially entering one polar cap or the other and thus explaining the 
north-south asymmetry of the polar rain. We model this anti-sunward moving 
strahl as a third bi-maxwellian 

f s = V<» 3/2 % » s 2) e *P<-v{;/u s ^ - v A V) <v„>0> (3) 

i " i 

where <j_ is considerably greater than u . 
s„ s 

To map an initial distribution function at point 1 (e.g. the solar wind) 
along a field line to point 2 (e.g. the tail lobe), we assume conservation of 
energy 

(1/2)mv^ + (1/2 )mv^ - e« 1 = (1/2)mv 2 + (1/2)mv 2 - e$ 2 (4) 

" X " JL 

and conservation of the first adiabatic invariant 


13 



( 5 ) 


kA /B 1 5 A /B 2- 

A 1 

Combining these equations we oan express the initial velocities at point 1 in 
terms of the final velocities at point 2 as 


A = <W V 2 

X X 

v« = v* + C1-« 1 /B 2 >v| - 2e/m (* 2 - 


( 6 ) 


A 

Invoking Liouville'3 theorem and assuming the field varies monotonically 
between points 1 and 2 to assure particle access to the stronger field, we 
know that f is conserved along the particle trajectory. Taking a velocity at 
point 2, we can find the equivalent velocity at point 1 from (6), evaluate f 
at point 1 from (1) and know from Liouville's theorem that it is the 
appropriate f for point 2. 

We have carried out this procedure and contours of f in the v„, v A plane 
are shown by the lighter curves in Figure 5. Figure 5a contours are 
calculated directly from equation 1 (neglecting f ) while those of 5b 

C 

represent the distribution mapped to a field which is 2.7 times the 
interplanetary field strength. We have assumed that $ 2 - =0 here and in 

subsequent discussions. In equation t a solar wind halo temperature and 

5 

density of 5x10 °K and ,1/cc was assumed along with a strahl temperature 
T„ a - 1.5x10^, T as s 2.6x10 J and n s = .086, These parameters were chosen to 
represent the data in Figure 2 and their selection will be discussed below. 

The oontours in Figure 5a are elongated along the v lt axis and are very similar 
to the strahl examples shown by Rosenbauer et al. [1976; 19773, Feldman et 
al., [1978] and Ogilvie and Scudder [1981] but are quite different from the 




more isotropic solar wind contours deduced from earlier wide angle solar wind 
instruments [e.g. Feldman et al,, 1975,1 which were widely known at the time of 
the Mizera and Fennell review [1970]. 

As the solar wind electrons move into the stronger field, conservation of 
y requires that their perpendicular speed increases in proportion to /D until 
all their energy is in v A , at which point they mirror and return to the solar 
wind. Only electrons from the shaded region in Figure 5a (those with pitch 
angles les3 than 37.5°) reach field strengths of 2. 7 B (Figure 5b) where they 
fill the entire v l|t v 1 plane. The contours in Figure 5b are noticeably less 
anisotropic even for this modest field increase. Those electrons that reach 
low polar cap altitudes experience a field increase of roughly 5000. They 
will be even more isotropic than in Figure 5b because they all come from a 
solid angle in the solar wind of less than 1°. There is not apt to be a very 
large anisotropy within suah a small solid angle and hence the low altitude 
distributions are expected to be much more isotropic. Note that the contours 
intersect the v„ axis at identical points in 5a and 5b indicating that the 
electrons moving precisely along the field (v^O) are unaffected by the 
changing field. 

Since most of the electrons passing a given point in the tail lobe will 
mirror at lower altitudes and move back out of the tail (excepting those 
within the very small and essentially unmeasureable loss cone of less than a 
few degrees) the fluxes should show little parallel-anti-parallel asymmetry 
relative to the field direction. This is what is observed by ISEE 1 in the 
tail lobes (e.g. Figure 2). 

Constant energy circles are overlayed on the distribution function 
contours of Figure 5a and 5b at the measured ISEE energies of 87 and 510 eV. 
Measurements made at various pitch angles along these iant energy circles 
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will show a decreasing f a3 the sampled pitch angle becomes oloser to 90°. 
Plots of log f versus pitch angle from the mapped model distribution are 
superimposed on the data of Figure 2 as smooth curves. Clearly both the 
absolute magnitudes and the pitch angle dependence of the model fit the data 
quite well. At high energies the model anisotropy is maximal due to the 
anisotropic strahl while at lower energies the distribution is more isotropic 
due to the increasing influence of the lower temperature, more nearly 
Isotropic halo. This lower anisotropy at lower energies is a very common 
feature in the tail lobes as it is in the solar wind [Rosenbauer et al. f 1976] 
and it thus supports the idea that the solar wind is the source of the polar 
rain. 

The technique used for determining the solar wind halo and strahl density 

and temperatures from tail-lobe data is illustrated with an example shown in 

Figure 6 and 7. The lower portion of the Figure 6 shows the log f versus 

pitch angle plots on March 10, 1978 in the format of Figure 2. At this time 

the experiment was in the higher energy mode extending to 7 keV but all 

energies above 1 keV were measuring background levels and those above 1.35 keV 

are not shown. At the top of Figure 6 we show plots of f versus v lf and v A 

which are the orthogonal auts through a three dimensional surface fit to the 

measurements [Fitzenreiter et al., 1 984 3. Using such fitted data is a slight 

improvement over the measured points since the measurements are made at pitch 

angles close to but not exactly at zero. In Figure 7 we plot the log of these 

2 

fitted f values versus v„. On this type of plot a maxwellian distribution 
appears as a straight line. At the highest energies we assume that only the 
strahl is important; we then can obtain the strahl parallel temperature 
{2x1 0^°K ) from the slope of this line and the coefficient of the strahl 
exponential from the intercept with the ordinate. For a maxwellian with a 
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single temperature this intercept would define n g , but with a two temperature 
bi-maxwellian it only defines a relation between n and T , and further 
information from other pitoh angles must be utilized to specify n g or T gA , In 
this case we find n g a .012 and T gA a 1.3x10 5o K. In principle the same 
technique can be used at lower energies where the halo dominates, but 
sometimes, as in Figure 7, there are not enough energy channels to define this 
slope. However, the pitch angle dependence of the data together with the 
limited energy interval may be combined to determine appropriate halo 
parameters (T h „ = 5.5x10 5 °K, T^ A = 1.5x10 5 °, and n h = .1/cc), A low energy 
slope determined in this way is shown on Figure 7. Winningham and Heikkila’s 
[1974] polar rain spectrum is superimposed for comparison and it shows a very 
similar temperature. The values from the model distribution mapped to the 
tail are also superimposed on Figure 6a and the very good agreement is 
apparent. These temperatures only characterize the more field aligned portion 
of the interplanetary distribution (equivalent to the shaded portion in Figure 
5a) and hence may not give a value of T A equivalent to that determined from 
the entire interplanetary distribution. The T )t halo values are, however, 
typical of the solar wind [Feldman et al., 1975]. This agreement between 
model parameters determined in the tail and typical 3olar wind parameters 
strongly supports the concept of solar wind entry. 

About four hours before the data in Figure 6a was taken, ISEE had 
observed the tail lobe at lower altitudes where the field strength was about 
100 nT or twice the value in 5a. This data is shown in Figure 6b. The curves 
from the model are again superimposed using the same solar wind parameters 
determined from 6a, correcting only for the greater lobe magnetic field 
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strength. The good agreement can be attributed to the validity of the mapping 
procedure and the stability of the souroe conditions over the intervening 
time. In neither case were field aligned potentials invoked. 

We note that the strahl temperatures determined in the above manner are 
typically 2x10^°K. This corresponds very well to the canonical value of the 
solar corona as would be suggested by an exospheric treatment [Schulz and 
Eviatar, 1972] and thus tends to further substantiate our proposed origin of 
these electrons observed in the tail lobes. We discuss this point more fully 
in section 6. 

5. Energetic Polar Cap Electrons 

One of the most interesting aspeots of polar cap electrons are the les3 
frequent occurrences of spatially uniform, higher intensity precipitation 3een 
at energies up to about 10 keV. As was mentioned above, this phenomena has 
generally been considered a3 a separate class of event which requires a field 
aligned magnetospheric acceleration mechanism. However, in view of our 
proposed association between the solar wind strahl and the polar rain, it is 
of interest to reexamine this problem and ask whether the characteristic 
variations in the strahl which have been reported [Rosenbauer et al. 1976; 
1977; Feldman et al., 1978; Pilipp et al., 1981; Pilipp, private 
communication, 1984] might explain these low altitude observations and 
eliminate the need for magnetospheric acceleration. 

Existing information on high intensity polar cap electrons ha3 come from 
4 events of roughly a day or two duration centered on December 20, 1971, March 
10, 1972 [Foster and Burrows, 1976], September 18 1974, and October 15 1974 
[Meng and Kroehl, 19771. All four events show intense precipitation 
preferentially on the northern polar cap and all events for which 
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interplanetary data is available (the latter three), are associated with away 
sectors. Also all the events are associated with magnetic atoms; Foster and 
Burrows [1976] note that their events are associated with post storm quieting 
and to some extent the same appears to be true with the events of Meng and 
Kroehl [1977]. 

Available solar wind data for these events is presented in Figure 8. 

Solar wind speed, density and proton temperature [King, 1979; 1983] are 
plotted versus time and shading indicates the approximate time and duration of 
the polar cap events. The March 10, 1972 event (Figure 8a) occurs as the 
speed decreases a few days after the arrival of a high speed stream on March 6 
which wa3 probably caused by a transient event on the sun. The polar cap 
event is closely associated with unusually low density solar wind of les3 than 
1/cc which caused the earth’s bow shock to move out beyond the IMP-6 
spacecraft which was located more than 25 upstream of the earth [Fairfield 
and Feldman, 1975). The September 1 97^ event of Meng and Kroehl (1977) is 
characterized by large polar cap ratios on September 16 which then decreased 
before increasing to even higher values on September 18. These two intervals 
are identified separately in Figure 8b. The first interval closely follows 
the arrival of a corotating solar stream and the second more intense event is 
associated with the latter stages of a very low density region as was the case 
on March 10, 1972. The October 15 event (Figure 8c) oocurs one solar rotation 
later and is associated with a corotating stream. The density averages 6/cc 
during this interval. 

ISEE 1 spent only 89 hours (Figure 4) in the tail lobe in 1978 when 
conditions were appropriate for the solar wind strahl entry to the northern 
hemisphere. Consequently it is not surprising that ISEE did not see any polar 
rain events as intense as those selected from the larger DMSP and ISIS data 
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sets. The example shown In Figures 6 and 7 is In fact a relatively intense 
event for ISEE and meaningful counting rates were seen as high as 761 eV. The 
interplanetary density was below Vca during most of this day. 

For these few examples then, it appears that the solar wind tends to have 
unusually low densities at times of more intense polar rain. Note that this 
observation is oontrary to what one might naively expect - that increased 
solar wind density would produce increased polar cap precipitation. We will 
return to this point after reviewing ideas on the strahl formation in the next 
section* 

Other important information about intense polar rain can be deduced from 
the shape of the distribution function. We have scaled points from Figure 3 
of Foster and Burrows [1976] and Figures 5 and 7 of Meng and Kroehl [1977] and 
these points, converted to distribution function f, are plotted versus energy 
in Figure 9a. Again quasi -maxwellian regimes will appear as straight lines. 
For reference we reproduce the two fitted distributions from our relatively 
intense ISEE example of Figure 6-7 where counting rates fall below the 
instrumental background above 1 keV. The shape of the distribution functions 
for intense polar rain events shows a general resemblance to the ISEE data 
that is characteristic of more normal polar rain. The Foster and Burrows 
data, which extend to low energies, have a low temperature component which is 
more intense than normal polar rain but has a '‘temperature" 

dlogf/dE =-1 /k"T" (7) 

similar to low energy ISEE polar rain and the solar wind halo. Between 200 eV 
and 700 eV all the intense events have high intensities but they have strahl- 
like "temperatures" of roughly 2 x 10 6 °K which are comparable to the ISEE 
data. The additional feature revealed by the intense events is a high energy 
tail above approximately 1-2 keV. We note that the intense polar rain spectra 
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are consistent with being proportional to the more typical ISEE polar rain, if 
we assume that the high energy tails for the intense events are due to 
intensity changes which raise the >2 keV tail above particle detector 
threshold. In the next section we review current ideas on the escape of 
electrons from the sun before concluding that intense polar rain data are 
compatible with a 3olar origin hypothesis, 

6, Electron Esoape from the Sun 

The current theoretical understanding of the formation of the 3 olar wind 
strahl suggests that it is shaped en route from the sun by the competition 
between two forces. On the one hand coulomb scattering broadens the 
distribution while on the other hand the Lorentz force, which in the absence 
of collisions is normally discussed via conservation of the first adiabatic 
invar lent jj, tends to induce a field aligned collimation [cf Jockers, 1970], 
The parameters that determine the importance of these forces are the flux tube 
profiles of density, temperature, and magnetic field strength. The coulomb 
scattering frequency v£v 0 (n/n o ) (T o /T) J ' <E 0 /Er is strongly energy dependent 
whereas jj conservation implies that the logarithmic rate of change of pitch 
angle is essentially independent of energy. Together these facts imply that 
there is significant energy dependence in the relative importance of the 
forces that shape f near 0° pitch angle. At sufficiently high energies the 
magnetic focusing makes f very nearly field aligned because coulomb collisions 
cannot broaden it. At low solar wind "core" energies the coulomb collisions 
dominate [Scudder and Olbert, 1979a] and the distributions are nearly 
isotropic. At intermediate energies there is a gradual transition and closer 
competition among the two effects. 


21 




i x i..su . _ . , j 

j 

! 

i 

f 

! 

I The approximately field aligned (v^O) "strahl" regime of the halo is 

| most easily thought of as the least coulomb screened remnant of the inner 

| coronal regions as illustrated by Soudder and Olbert [1979a, eqtn 113 where 

] 

j (in the high energy limit) 

j f # (E^, M = °) = W ^ exp(-S(E^, E q , n, T)) + low energy correction (8) 

Here f # is the distribution function at earth, while f 0 is the corresponding 

quantity in a region of rapidly decreasing density at 1.03 R . The appropriate 

s . 

! form of f is not known a priori, but it surely depends on the density and 

scale heights in this vicinity. At times when coronal densities attain higher 
or lower values, the function f Q will in general change. In equation 8, S is 
an energy and path dependent counter of the number of collisions a fictitious 

i 

! "collisionless" particle would have encountered getting to the observer: S = 

/v(E (t) ,t)dt, where v is the collision frequency at time t along the particle 
trajectory. Note that the energy E(t) of the electrons leaving coronal 
altitudes is reduced by climbing up the heliospheric potential barrier. This 
raises the local collision frequency and makes S a spatial integral over the 
local collisional frequency and the potential between the sun and earth. The 
coulomb terms become increasingly important as the local kinetio energy is 
reduced and or as the effective path length is increased due to the pitch 
angle increasing. 

The energy dependent counter S determines the opacity of the coronal 
distribution function f 0 to the local observer. This opacity is determined by 
the radial profiles of the electron density, temperature, and heliospheric 
potential. In the high energy regime, where (8) is approximately valid, S’ 
affects f in a simple way given by 

log 10 f # = lo 8^ 0 f o ( n * E > + ln(correction) , (9) 
where S’ = Sxlog 10 e, is the energy dependent number of decades f Q is reduced 
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to produce f . Figure 9b illustrates S* (n,E ) which has been calculated for 
two different coronal density profiles corresponding to a polar coronal hole 
[Munro and Jack3on, 19773 and an equatorial region. (Parameters for these 
profiles are from Olbert [private communication, 1978] and correspond to 1 AU 
densities of *1. 8/cc and 15. 4/cc. The electron temperature ha3 been assumed to 
vary as n 0 which implies that the heliospherio potential energy drop between 
point r and infinity is 7kT(r) which distills the results of Sittler and 
Scudder [1980].) The quantity S* is large at energies below the heliospheric 
potential energy (=750 eV) where scattering dominates but as E increases, S 1 
rapidly beoomes very small for both profiles. Although S’ never attains a 
value of zero indicating a complete absence of screening, we adopt a value of 

S' = 0.5 as an effective lower bound on a high energy region where screening 

* 

is not important. Above this energy, E , the spectral shape at 1 AU is little 

modified from f . Figure 9b shows that E is higher for the higher density 
o v 

profile. Because of the exponential screening factor this leads to the 
important expectation that higher (lower) column densities should strongly 
depress (enhance) strahl phase densities. 

The portion of f that reaches 1 AU is shown explicitly in Figure 9c. 

Here we have taken the Foster and Burrows polar cap spectrum from Figure 9a 
which was obtained at a time of very low solar wind density. The top curve is 
3pline approximation to the essential features of this spectrum. The heavily 
shaded region plus the diagonally shaded region of this figure represents that 
portion of f calculated as surviving at 1 AU if the heavy curve were assumed 
to be the actual coronal boundary condition and realistic low density 
conditions prevail. Above 200 eV the observed and primary spectrum are nearly 
identical and in the usual polar rain regime 200-500 eV range the slope is 
essentially that of f in the corona. The diagonally shaded region indicates 
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what would remain from the same solar spectrum if the column density were 
high. This diagonally shaded region is strongly depressed relative to the low 
density case indicating that the observed distribution in the usual polar rain 
regime has a slope that is severely modified. In fact the actual slope in 
this region (not shown) will be that of the nearby (cooler) halo which will 
fill in the low energy void (i.e., the low energy terms neglected in equation 
8). Notice that the slopes of the data ourve3 (Figure 9a) in the 200-500 eV 

g 

range are most nearly coronal (2x10 ) for the low density and or high speed 
examples in Figure 9a and cooler for higher densities. 

The above theoretical concepts are in agreement with much of what i3 
known about the strahl, Experiments (generally at E<1 keV) have shown that the 
strahl i3 most prominent (intense and narrow) in high fpeed streams 
[Rosenbauer et al., 1976; Pilipp et al., 1981] where the density is on average 
low Feldman et al., 1978], At 1 AU the solar wind number flux, nV, tends to 
be constant on average (Bridge, 1976), although this relation need not 
strictly hold in individual cases. .Tri comparisons with the data that follow, 
we will use thu locally measured density as the observable most closely 
related to the integrated column density. At the same time we will remember 
that local density enhancements may oocur and that high speeds may be an 
additional indicator of low column densities. 

Little is known about the strahl in the solar wind above 10 keV since its 
narrowness makes it is very difficult to detect. Strahl '•temperatures*' (in 

the sense of equation 7) at solar wind electrostatic analyzer energies (E =300 

5 5 

eV) have been reported between 7x10' [Ogilvie and Scudder, 1981] and 8.9x10 

[Feldman et al., 1978] and a3 high as 2x10^ [Rosenbauer et al., 1976], The 
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latter authors noted that no single maxwellian character ized the entire energy 
range in the strahl direction with the highest T estimates (in the sense of 7) 
associated with the highest halo energies (<1keV). 

The labels of the intense polar rain speotra in Figure 9a have boon 

arranged in order of increasing, 1 AU solar wind densities. At one extreme in 

Figure 9a the two most intense polar rain speotra in the 200-500 eV range 

correspond to an unusually low solar wind densities of 0.6 and 0.9/oc. The 

fact that no solar wind antecedent could be found for the Foster and Burrows 
partioles by wide aperture solar wind electron detectors is consistent with 
the theoretical prediction that maximum field oollimation (producing maximum 
difficulty of detection) would be expeoted at this time. At the other 
extreme, the Fennell et al„ sample which revealed similar polar aap and solar 
wind intensities occurred when the solar wind density was upwards of 20/co. 
Under these conditions the field aligned component is theoretically expected 
to be significantly depressed and broadened. Yhe other high intensity cases 
abstracted from those of Heng and Kroehl either occur at similarly low 
densities or aooompany relatively high speed solar wind >550 km/s. Together 
these examples reveal an association of high intensity polar rain and 
extremely low density and or high speed solar wind at 1 AU. 

Although column density variations can explain the strahl intensity 
changes at intermediate energies (E <200-500 eV) where electrons are scattered 
away from field alignment to a greater or lesser degree, at higher energies 
field alignment always prevails and hence the variations of S’ cannot explain 
the higher energy intensity changes. At these higher energies (where S* <0.5) 
the 1 AU observer sees an essentially unscreened f o ; the question thus beoomes 
how doe3 f^ vary near 1,03 R at these times. Since high intensities appear to 
be associated with low densities and high speed streams we consider solar 
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conditions at these times. These wind conditions are widely held to result 
from extremely low density and or rapidly diverging magnetio topologies over 
coronal holes [Holzer, 1977} Levine et al, 19773. These regions also have 
stronger density gradients than regions where slower flows appear to originate 
[Munro and Jackson, 19773. 

It is indeed reasonable that f depends on the absolute coronal densities 
and scale heights at 1.03 R Q . In particular, in the presence of strong 
gradients in a fully ionized plasma at the surface of a star it is very 
difficult to guarantee a priori that a Maxwell Boltzman distribution would 
adequately describe the highest looal kinetic energies. In fact it is very 
difficult to surpre33 the non-thermal tails [Soudder and Olbert 1979a, b, 

19833, especially in the steep ooronal hole density profiles which also have 
low absolute density. These tails usually appear at approximately 7 times the 
local most probable energy partiole [Soudder and Olbert, 1 97 9* a b3 which at 
coronal temperature of kT =170-200 eV implies significant, flatter non-thermal 
tails above break energies of 1190-1400 eV at the 3un. At the orbit of earth 
these electrons have given up kinetic energy to the heliospheric potential 
(500 to 1000 eV), Modest net energ*; rh* fts are anticipated aoross the bow 
3hock-magnetopause so estimates of net energy shifts for the breaks in the 
injection spectrum at earth are between 425 and 1000 eV. This regime is 
compatible with Figure 9a which shows a transition from (-kT ) _1 =dlnf/dE =1.5- 
2x10^°K in the regime 200-300 eV with transitions to flatter tails occurring 
at energy from j* 500 eV for Foster and Burrows to 1300 eV for Meng and Kroehl. 

An exciting aspect of polar solar measurements is that the magnetic 
defocusing of the strahl over a oollisionless traversal of the shock, sheath 
and and lobe, allows a much better view of theory alonB B beaause it takes 
the narrow bundle of phase space and disperses it over a much wider solid 
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angle 30 that f Qb3 (<v A >=0) in the magnetotail or even better at lower 
altitudes is essentially theory ( v j = 0), a quantity most direotly related by 
exospheric theory to the coronal boundary conditions f 0 [Schulz and Eviatar 
1972], Recall that in these exospheric theories, there is no guarantee that 
the boundary condition of f Q equal to a Maxwellian is physically warranted. 
Polar cap measurements of the high energy form of f can help establish f Q and 
contribute to the longstanding issue of how to energetically support the high 
speed solar wind state seen at 1 All. 

The detection of the supra-coronal thermal spectrum has profound 

implications for indirectly determining the coronal temperatures from the 

interpretation of optical emissions from the coronal regions and the freezing 

in of coronal charge states sampled in the solar wind [Owocki and Scudder, 

1983} Owocki, 1982]. The existence of strongly suprathermal distributions at 

the base of the solar wind expansion has also been previously suggested by 

Scudder and Olbert [1983] and Olbert [ 1 962 ; 1 983 J as being required for a 

conduction supported high speed solar wind. Finally it is of interest to 

estimate the implications of the polar cap observations of f to determine the 

implied temperature of the coronal electrons (in the sense of P/nk) and 

contrast it with the upper bounds of T=2x10^°K suggested by Leer and Holzer 

[19793. Assume a two maxwellian approximation for f # above E=100 eV. Let T co 

be the effective "temperature" for 100 eV< E # <400 eV and T^^ the 

corresponding quantity for the flat tail regime.. Let their corresponding 

densities be n Q ° and n h ®. Let the potential energy difference to r=1.03 R° be 

given by ea#. Using the values fit to Foster and Burrows in Figure 9 yields 

n h /n =0.05; T ®=10T assuming e.A»=750 eV implies that T =1.009 T or that 
n c n c 0 co 

monitoring T ue essentially monitor T . Thus the presence of the 
suprathermal tails, a few percent by number irj the corona with a high 
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"temperature”, does not contradict the ooronal bounds of and are consistent 
with arguments of Leer and Holzer [19793. 

7. Summary and Discussion 

In this paper we have presented new evidence that polar ra*n electrons 
observed in the magnetotail are highly anisotropic in the field aligned sense. 

We have used this stringent new observational constraint to test the solar 
wind entry hypothesis. We find support for this idea based on the 
quantitative agreement between lobe measurements and typical solar wind 
distribution functions mapped to the tail lobes. Furthermore we find that it 
is the solar wind strahl portion of the 3uprathermal electrons which is 
responsible for both the tail lobe anisotropy and electron precipitation over 
a preferred polar cap. Alternative explanations fail to provide such ; 

agreement; the more traditional and more isotropic core and halo solar wind 
electron distributions known at the time of Fennell et al. [19753 would become 
even more isotropic when they enter the tail and hence they fail to explain f 

the observed lobe anisotropies, ii 

Although acceleration by a field aligned electric field might seem like a j 

method for introducing such an anisotropy, this does not appear to be a valid | 

alternative. A field aligned electric field along an intensifying magnetic 

ie 

field does not preferentially accelerate field aligned particles within a hot jl 

ii 

population. Although field aligned particles are accelerated, those with [j 

slightly larger pitch angle are also accelerated and isotropy tends to be p 

H 

preserved. Only if one accelerated a very cold distribution with low v, [j 

ii 

spread can one preferentially produce field aligned electrons and avoid the 

2 ! 

accompanying less field aligned electrons as B increases and v A /B is i 
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conserved. If the solar wind core electrons were accelerated to give the 
observed anisotropy, the phase space density and temperatures would not matoh 
the tail observations, 

7.1 Magnetospheric entry 

We summarize our views of electron entry to the magnetosphere using 
Figure 10 which schematically illustrates electron distribution funation 
contours relative to the magnetio field in various spatial regions. All 
distributions are supported by solar wind and magnetosphere measurements 
described above or magnetosheath measurements discussed below. 

Distribution 1 in undisturbed interplanetary space is the halo-strahl 
distribution with the asymmetry provided by the almost oollisionless strahl 
electrons streaming out from the sun. Note that the halo-strahl distinction 
is a convenience for modeling and describing a distribution funation that 
gradually changes from highly anisotropic to isotropio with decreasing energy; 
in reality the low energy isotropy is probably due to the highly-scattered, 
outward-moving coronal electrons that would be highly anisotropic in the 
absence of collisions. Hie shaded area along the field direction 
schematically illustrates those electrons destined to precipitate on the polar 
cap if the field line connects to the tail. In distribution 1 the actual 
angular width of the region is of the order of 1°. Distribution 2 is the 
field aligned population measured by ISEE 1 on tail lobe field lines that 
connect to field lines leading back to the sun. The distribution is field 
aligned due to strahl electrons freely entering the tail but it is less 
anisotropic than in the solar wind due to the effect of the electrons moving 
into the stronger lobe field. The distribution is bi-directional because most 
of the electrons passing ISEE mirror at lower altitudes and return to the 
spacecraft to produce the nearly symmetric distribution. Only the most field 
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aligned electrons within about 2° of t) precipitate and produae the polar rain. 
The upcoming loss cone 3hould be devoid of solar wind electrons because these 
particles have precipitated. Thi3 region may be populated by ionospherio 
sources, but it i3 too narrow to be measured by ISEE. At low altitudes over 
the polar cap the distribution is much more symmetric and the loss cone fills 
most of the down coming hemisphere (distribution 3). 

In the other tail lobe where the reconnected field lines lead to 
interplanetary lines that do not lead directly back to the sun, only the more 
isotropic sunward half of the interplanetary distribution (the halo) enters 
the tail. This distribution mirrors and produces distribution 4 which is 
isotropic except for the empty loss cone. ISEE only measures these isotropic 
electrons at low energies in the "unpreferred" polar cap. The quasi-isotropic 
halo is seen over both polar cap3 but over the preferred cap it is 
supplemented by the anisotropic electrons of the strahl. The presence of the 
energetic strahl electrons over one pole but not the other leads to large 
north-south asymmetries at high energies. The more isotropio low energy 
distributions over both polar caps leads to decreasing north-south asymmetries 
at these lower energies. Changing solar-interplanetary conditions will raise 
or lower the energy where the north-south asymmetries become large. 

In the magnetosheath an anisotropic distribution 5 would be expected on 
field lines that simply continue on downstream of the earth. This 
distribution would be reversably broadened but its field aligned asymmetry 
preserved. Where solar particles traverse the thinnest perpendicular portion 
of the subsolar bow shock they remain magnetized for gyroradii r g <<15 km which 
is the thinnest scale of the perpendicular shock [Scudder et al.1984]. This 
implies that electrons with energy below ^500 eV only gain energy from the 
deHoffman-Teller frame electrical potential which is =30-80 V [of Goodrich and 
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Scudder, 1984]. Higher energy electrons with larger gyroradii could gain 
energy of up to 200-300 eV across a perpendicular shook, but less across a 
qua3i-parallel shock. In reality, most polar rain electrons will approach the 
more distant magnetopause having traversed a weaker distant bow shock where 
they will have gained much less energy than the worst case numbers above. 
Distribution 6 or 7 ought to be seer, on magnetosheath field lines that connect 
to the tail.. Which of these distributions is seen will depend on which polar 
cap the distribution is mirroring in. 

The anti-sunward moving half of the anisotropic magnetosheath 
distribution (5 or 6) has probably been seen by the Apollo plasma detectors on 
the surface of the moon [Reiff and Reesoner, 1975]. They frequently detected 
two maxwellian components in the magnetosheath with the hotter component 
having temperatures and densities of 2.5x10^ and .05/cc which are similar to 
the then unknown strahl. This component was anisotropic (different 
intensities in two detectors looking in different directions) and time 
variable. Although these authors suggested bow shock electron acceleration as 
the cause, in retrospect it seems quite possible that a varying magnetic field 
was moving the anisotropic strahl in and out of the fields of view of their 
detectors. 

Since we advocate the entry of solar wind electrons to the tail lobe, it 
might be thought that the low energy, high density electrons should also enter 
the magnetotail. Such entry would be contrary to longstanding evidence for 
low density lobes, an observation confirmed indirectly on ISEE by the 
existence of a spacecraft plasma sheath that would not occur in the presence 
of a higher density plasma. We suggest that the lack of intense low energy 
electrons involves the behavior of protons. The low energy electrons, having 
much higher speeds and greater mobility than protons, could easily follow the 
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field lines into the magnetosphere, but to do so would produce field aligned 
currents and an unacceptable charge imbalance. We propose that a field 
aligned thermal potential is created near the tail magnetopause to prevent 
suoh low energy electron entry. This potential drop will ocour across the 
strong density gradient in the reconnected flux tube discussed by Sibeck and 
Siscoe, [1984], This potential will adjust itself to allow the appropriate 
number of electrons to enter consistent with the rotational discontinuity mass 
flux requirements. 

We emphasize that this magnetopause potential decelerates the entering 
electrons whereas previous work (e.g., Foster and Burrows, 1977; Gussenhoven 
et al,, 1984) has concentrated on accelerating tail lobe electrons. In our 
picture the bow shock accelerates the electrons while the magnetopause 
decelerates them to a somewhat greater extent. This situation leads to tail 
lobe energies slightly reduced relative to the solar wind that are still 
consistent with our earlier assumption of 41=0 made while comparing ISEE 
measurements with typical solar wind strahls. 

Furthermore we suggest that the magnitude of this magnetopause potential 
will be related to the location where the field line enters the magnetopause. 
On magneto spheric cusp field lines the potential as a fraction of the 
magnetosheath electron temperature will be small and large fluxes of hot, 
magnetosheath electrons and their accompanying protons will enter as is 
observed [e.g., Winningham, 1972], On field lines poleward of the cusp the 
ions cool faster than the electrons, and simultaneously the magnetosheath 
density is reduced as the flow increases. Both these changes imply that 
increased decelerating potentials are required to permit ever smaller fluxes 
of electrons (and ions) to move down lobe field lines. A calculation of this 
potential profile demonstrates that the equilibrium electron flux (mean 
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energy) is most intense (lowest) near noon and is more than one order of 
magnitude larger (factor of 3 lower) than the fluxes (mean energy) near 
midnight, Thi3 noon-b, '(-midnight variation in intensity is anti-aorrelated 
with energy and is the main experimental finding in the reoent paper of 
Gussenhoven et al. [ 1 984 3. These authors suggested that internal 
magnetospherio processes were the probable cause of this morphology. Our 
suggestion is that the polar rain electrons come directly from the coronal 
solar wind with magnetopause boundary condition systematics modulating their 
mean intensity and energy. Energization need not occur anywhere within the 
magnetosphere. 

7.2 Solar Source 

Once the strahl wa3 established a3 the source of polar rain electrons, we 
asked the further question of how the strahl is produced by the sun. Current 
thinking is that the strahl is due to the nearly collisionless escape of 
electrons from the corona, with the collimation along the field direction 
produced by adiabatic focusing in the decreasing solar-interplanetary magnetic 
field. An energy dependent collision frequency dictates that different 
energies will arrive at 1 AU with different anisotropies. 

Lower energy electrons should undergo frequent collisions and be severely 

scattered. In this manner the quasi-isotropic halo will be produced from 

; 

electrons that would otherwise have been field aligned. At slightly higher 
energies there ought to be a closer competition between scattering and the 
tendency toward field alignment. Lower solar wind column densities will tend 
to make field alignment prevail and produce higher field aligned intensities 
that extend to lower energies. Effective "temperatures" in this intermediate 
energy range should correspond to coronal values. For higher solar wind 


densities, scattering will prevail, reducing field aligned intensities, and 
shaping a more isotropic, cooler, halo-like distribution. 

At still higher energies, electrons with lower collision frequencies 
should make the transit to 1 AU with little scattering regardless of column 
density. They will become very field aligned at 1 AU and preserve the true 
characteristics of the high energy spectrum at the solar source. At these 
higher energies we must look to variations at the solar source to produae 
variations at 1 AU. We noted that times when the low column densities will 
produce enhanced strahl formation tend to be associated with high speed 
streams and low densities. These conditions are thought to originate from the 
central regions of coronal holes where the density is low and the scale 
heights of density and magnetic field are shorter than in the source regions 
for the slow, dense solar wind. It is, therefore, more difficult to suppress 
suprathermal tail formation in these regions than elsewhere. Accordingly, it 
is expected that f Q above source regions for fast, sparse solar wind flows 
should be preferentially more non thermal and intensities at higher energies 
should be enhanced. This implies that an observer at 1 AU who is essentially 
unscreened from the corona at sufficiently low energies may see signatures of 
the statistical equilibrium in the corona, including the production of a 
non-thermal distribution with a "temperature" which is nearly 2x10^°K. 

7.3 Energetic Polar Cap Precipitation 

The above concepts are in fact supported by observations. First of all 
an enhanced strahl is associated with high speed streams which in turn are 
known to be associated with low densities. We also noted the similarity of 
the f versus E plots in Figure 9a for the intense events and the more normal 
polar rain. Of particular note was the common =2x10 6 °K coronal-like 
"temperature" (i.e. slope) in the energy range =200-700 eV. In fact the 
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primary difference between the normal and intense events seems to be 
proportionately higher fluxes and the presenoe of a high energy tail above 
about 1 keV during the intense events. This high energy tall could well be 
present during the more normal events, but at these lower intensity levels it 
would be below instrumental background. Finally we noted that the intense 
events seem to be associated with intervals of an usually low 3 olar wind 
density which, as was pointed out above, is the 1 AU condition for pronounaed 
strahls . 

Arguments against a solar origin given by Foster and Burrows [1976] were 
(1) the non detection of high fluxes of interplanetary source electrons at the 
time of intense polar cap precipitation events, and (2) the similarity of the 
spectrum of polar cap and plasma sheet electrons whioh they felt implied a 
common source. With modern information on the observed narrowness of the 
strahl and the expectation that 3till higher energy electrons would be even 
more collisionless and field aligned, it now seems quite likely these 
electrons would have eluded detection by the wide angle detectors available at 
that time. As for point 2, we agree that the two spectra are similar, but the 
solar wind/corona could equally well be the common source for these two 
regions on March 10 1972; the electrons could have E x B drifted onto closed 
plasma 3heet field lines or they could have become trapped on field lines that 
became closed due to tail reconnection. 

7.4 Conclusions 

Based on the quantitative mapping of modern 3olar wind measurements to 
the tail lobe we conclude that normal polar rain is of solar origin. We 
explain existing measurements with a solar wind entry model that doe 3 not 
require accelerating field aligned potentials. A more rigorous solution to 
the complete polar electrodynmics problem that considered atmospheric 
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photoeleotron3 and polar wind populations might well require weak field 
aligned potentials, but their effect on the polar rain at energies £100 eV is 
probably minimal. Unusually intense polar rain seems to appear at times when 
the strahl is expeoted to be especially prominent and solar electrons are apt 
to be intense. Hence direct solar aoces3 seems to provide a simple, viable 
explanation for intense polar rain as an extreme of coronal boundary condition 
Variations. This straightforward explanation is to be contrasted with 
magnetosphenc acceleration theory [Foster and Burrows, 1977] whioh requires a 
number of untested assumptions. Although we prefer solar origin as an 
explanation for intense polar rain, there remains an unproven possibility that 
spatially confined processes that undoubtedly work in the magnetosphere to 
produce polar oap arcs might operate on broader soale3 and produce the 
non-uniform intense precipitation. Even if the latter case is advanced, we 
feel that solar wind boundary conditions as outlined here remain important in 
determining (1) the original particle source and (2) the extent of the 
required particle acceleration. 

8. Future Work 

The access of solar electrons to the earth's polar cap3 produces numerous 
opportunities for further studies in solar and magnetospheric physics. It is, 
however, important to understand the degree to which field aligned potentials 
affect measurements in different regions. Comparisons of polar cap and 
magnetotail fluxes along the same field lines should be able to determine the 
significance of acceleration within the magnetosphere. Comparisons with 
interplanetary and magnetosheath data remain highly desirable in determining 
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the effects of bow shook and magnetopause potentials, but will have to await 
future instruments that are able to measure highly collimated electron 
distributions. 

Polar rain measurements are of great potential value for solar physics. 
The problems of directly measuring interplanetary fluxes within one or two 
degrees of the magnetic field direction can be attacked by using the earth’s 
gradually increasing field to broaden the field aligned beams to measurable 
widths. For instance, Foster and Burrows [1977, figure 73 measured a 4 keV 
intensity ratio of 1.9 between 0° and 60° at 1400 km altitude over the polar 
cap. Adiabatic mapping implies the same variation is present between 0° snd 
1° in the solar wind. This is an even more peaked distribution than our 
strahl model but it cannot be directly tested with existing solar wind 
measurements because observations in weak fields cannot begin to resolve such 
narrow pitch angle distributions. Such structure is, however, predicted by 
Olbert [19833 even when coulomb effects are considered. Interpreting such low 
altitude pitch angle data in this manner could yield much information on the 
solar origin and interplanetary propagation of such particles and contribute 
to the ongoing study of the problems of accelerating the high speed solar wind 
[Leer et al., 19823. 

That part of the polar rain distribution that is determined by collisions 
in the lower corona contains information on coronal temperature. This is an 
important, but poorly known solar parameter and polar rain measurements could 
potentially be the best way of directly determining this quantity and the 
state of statistical equilibrium that exists at the coronal base. 

In the magnetosphere the polar rain may well provide a clear and 
definitive measure of open field lines. Energetic solar electron data have 
been interpreted in this manner for many years [e.g. Paulikas, 1974] but such 
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energetic particles achieve measurable flux levels only when they are produoed 
by large solar flares. With more continuous lower energy data one can hope to 
monitor the polar cap size end shape and study its variability and its 
relationship to regions such as boundary layers. 
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Figure Captions 

Figure 1 Schematic Illustration of how interplanetary magnetio field lines 

conneot to those of the earth. In a ''toward 1 ' interplanetary seotor 
(top) the interplanetary field lines that lead baok to the sun 
connect to the southern tail lobe and polar cap whereas in an "away" 
seotor (bottom) the solar connection is with the northern hemisphere. 
Open arrows indicate the motion of electrons coming from the sun. 

Figure 2 ISEE~1 pitch angle distributions illustrating the nature of 

anisotropic tail lobe fluxes. The log of the distribution function 
is plotted versus pitch angle with every traoe offset by 1,25 deoades 
to prevent overcrowding. The instrumental one count level is 
indicated by a dashed line. The effective background level is 
slightly above the one count level especially in the sunward looking 
hemisphere between 90° and 150°, Meaningful counts are concentrated 
near the field direction at energies above 180 eV. Smooth curves 
represent a model solar wind distribution function mapped to the 
magneto tail. 

Figure 3 A six hour interval of data where anisotropic lobe electrons cause 
f„ (lines) to be above f a (dots) at energies below 1 keV. Log f is 
plotted for each energy with adjacent energies offset by 1 decade to 
prevent overcrowding. Enhanced energetic electron fluxe3 after 0700 
indicate entry to the plasma sheet. 

Figure 4 Histograms of f„/f A whioh indicate how larger 180 eV anisotropies in 
the northern tail lobe are associated with away interplanetary 
sectors. These average f ratios are reliable indicators of the 
anisotropy even though they drastically underestimate the magnitude 
of the anisotropy (see text). 
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Figure 5 Contours of a model electron distribution function (light lines) 


constructed to represent data of Figure 2 , Figure 5a represents the 
distribution in the solar wind and 5b illustrates the 3ame 
distribution when mapped to the stronger fields of the magnetotail, 

A decreasing degree of anisotropy can be seen in going from the solar 
wind to the magnetotail; a detector sampling along the 5)0 eV oircle3 
in the solar wind would sse f decrease by more than 6 orders of 
magnitude whereas in the magneto tail the decrease would be less than 
3 orders of magnitude. Only electrons within the shaded region of 
the solar wind distribution (5a) are able to reaoh a tail lobe field 
2.7 times larger (5b) where they form the entire distribution. 

Figure 6 Two ISEE-1 lobe pitch angle distribution' in the format of figure 2 
are shown in the lower portion of the figure. Data are from 
individual spacecraft spins and the low counting rates near 90 s are 
frequently at the one count level. The top portion of the figure 
shows plots of log f versus v where the 2 traces are parallel (v^sO) 
and perpendicular (v„=0) cuts through a three dimensional surface fit 
to the measurements. Dots indicate one count levels. Meaningful 
counts are again seen only along the field direction and only the f 
trace in the top portion is above the 1 count level. Smooth curves 
in the bottom portion represent the same solar wind model mapped to 
two different lobe field strengths.. 

Figure 7 The f data from tne top portion of Figure 6a are plotted on a semi- 

2 

log plot vs v so that straight lines indicate effective temperatures 
in different energy regimes. Tail lobe data above about 200 eV are 
generally consistent with a coronal like temperature whereas data 
below 200 eV are consistent with a solar wind halo temperature. 
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Figure 8 Solar wind plasma parameters for three published examples of intense 
polar rain. Low solar wind densities and or high speeds tend to be 
associated with the shaded regions of intense polar rain. 

Figure 9 (a) Illustrating various published examples of high intensity polar 
cap precipitaion along with the ISEE-1 fit from Figure 7 and an 
example of Fennell et al, taken during dense, slow solar wind 
conditions. In all cases the slopes in the energy ranges <200 eV and 
200-500 eV are approximately equivalent to solar wind halo and strahl 
temperatures respectively, (b) Indicating the energy dependence of 
the screening factor which determines how the 1 AU field aligned 
distribution function is related to the 3olar source distribution. 
Small values of the screening factor indicate that the source 
distribution is essentially unmodified at 1 AU whereas large values 
indicate that collisions and solar potential effects produce a vastly 
different distribution at 1 AU. High solar wind intensities (top 
curve) greatly enhance this modification, (c) Indicating a solar 
source distribution (top curve) and the portion of this distribution 
calculated as reaching 1 AU for the two different solar wind models 
used in 9b. A greater portion of the source distribution reaches 1 
AU when the solar wind density is low. An effective solar 
"temperature" can be deduced from polar cap measurements at these 
times . 

Figure 10 Schematic illustration of electron distribution function contours in 
various spatial regions in the vicinity of the earths magnetosphere. 
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Response to Referee 02 


Below are our responses to the five points raised by referee 02. The 
only other substantive changes involve references to a recent paper of 
Gussenhoven et al. and an explanation of their results via the magnetopause 
potential (p. 32), 

1. We have made numerous corrections to the equations. The two sets of 
dashed lines in Figure 9a have been clarified. Figure 10 has been 
corrected. 

2. We hope all the figures and references to them have been corrected. 

3. The high energy tail in Feldman’s data is at strahl energies. What was 
not appreciated at this time was the narrowness of the distribution. 

Indeed this anisotropy would produce north-south asymmetries and this had 
been suggested (as we stated in the introduction), but no quantitative 
comparisons with 30lar wind data had been attempted. We have inserted a 
sentence on p. 29 indicating that the halo-strahl distribution is 
convenient for modeling and describing the distribution but that there is 
no clear difference. 

4. No attempt has been made to search for the strahl in ISEE data in the 
magnetosheath or interplanetary medium. It would be difficult (though not 
impossible) to do, but this would be a whole new study. The effects of 
the bow shock are discussed on page 30-31. A sentence has been added 
emphasizing that most polar rain particles probably traverse a distant and 
weak bow shock that would have little effect on the particles. Figure 10 
is indeed based on the general knowledge of what the solar wind 
distribution function looks like. 

5. Apparently no ions have been detected along with electron precipitation. 

The plasma sheet appears to be clearly present at lower latitudes in the 
Foster and Burrows data and there appears to be little question that the 
observations are a real polar cap phenomena. The interplanetary B is 
only slightly positive during some of this time. z 
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